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The well-known coupled-mode theory of waveguides is extended to in- 
clude dielectric guides made of anisotropic materials. Exact coupled-wave 
equations for anisotropic dielectric waveguides are derived, and explicit 
expressions for the coupling coefficients are given. The coupling coefficients 
for isotropic waveguides are obtained as a special case. A simple approxi- 
mation for the coupling coefficients in the case of slight anisotropy and 
slight departure from an ideal waveguide is presented. 

I. INTRODUCTION 

The theory of dielectric optical waveguides deals with electro- 
magnetic wave propagation in optical fibers and in the waveguides 
used for integrated optics. Wave propagation in these structures is 
described in terms of normal modes. 1-3 However, normal modes pre- 
serve their identity only in perfect waveguides without irregularities 
of either the refractive index distributions or the waveguide geometry. 
Electromagnetic wave propagation in waveguides with any kind of 
irregularities must be described by means of coupled-mode theory. 3,4 
The electromagnetic waves in imperfect waveguides are expressed as 
superpositions of all the modes of a perfect waveguide. The mode 
amplitudes are coupled together by coupling parameters that depend 
on the nature of the waveguide imperfections. A description of wave 
propagation by means of coupled-mode theory allows calculation of 
radiation losses caused by intentional or unintentional fluctuations of 
the refractive index along the axis of the waveguide or by core-cladding 
boundary fluctuations. 23 Coupling among guided modes is used to 
design modulators or distributed feedback circuits for lasers or to 
effect improvements in the multimode dispersion properties of over- 
moded waveguides. The coupled-mode theory is well developed for 
waveguides that consist of isotropic dielectric materials. 3 ' 4 Some work 
has been done to extend this theory to waveguides consisting of aniso- 
tropic materials. 5-7 These waveguides are assuming increasing im- 
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portance in integrated optics as methods are being perfected for 
fabricating waveguides by diffusing different dopants (or outdiffusion 
of certain component atoms) into anisotropic crystals. 8-10 

This paper describes the derivation of coupled-wave equations for 
the modes of waveguides consisting of anisotropic materials. The 
coupled-wave theory is based on the definition of guided and radiation 
modes as solutions of Maxwell's equations for idealized structures. An 
orthogonality relation is derived that is needed to isolate individual 
terms in the infinite series expansion of the electromagnetic field. The 
principal result of this theory is the derivation of coupling coefficients 
that are important for solving coupled-mode problems. Readers not 
interested in the derivation should look at eqs. (46) and (48). Applica- 
tions of this theory are not presented here, since they will be the sub- 
ject of further publications. 

II. THE FIELD EQUATIONS FOR ANISOTROPIC MEDIA 

The derivation of coupled-wave equations for anisotropic dielectric 
waveguides follows closely the procedure used for deriving coupled- 
wave equations for isotropic waveguides. 3 The objective of coupled- 
wave theory is to construct solutions of Maxwell's equations for wave- 
guiding structures consisting of general refractive-index distributions. 

Anisotropic media are characterized by a dielectric tensor, 

(Cxx *xy *xz \ 
tyx *vv *yz )• (1) 

izx € 2 y €zz J 

We assume that the elements of this tensor are real quantities char- 
acteristic of lossless materials. It can be shown that conservation of 
energy requires that the dielectric tensor form a symmetric matrix 
so that the following relations hold : n 

txy = tyx; Cxz = £zx) tyz = €zy> \6) 

The magnetic properties of the medium are assumed to be the same 
as that of a vacuum so that we use the (isotropic) magnetic permea- 
bility constant no- Maxwell's equations for anisotropic media assume 
the form 

VXH = tco6-E (3) 

V X E = - to/xoH. (4) 

It was assumed that the electric field vector E and the magnetic field 
vector H have the time dependence, 

e iat . (5) 
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The tensor notation e-E may be expressed in component form as 

(6-E),- = euEj. (6) 

Summation over double indices is understood, and the subscripts i 
and j assume the values 1, 2, and 3 that represent the x, y, and z 
components of the vector E or tensor e. 

Derivation of coupled-wave equations for isotropic media is facil- 
itated by expressing the longitudinal components of E and H in terms 
of the transverse components. 3 This practice is preserved for our deriva- 
tion of coupled equations for anisotropic media. We single out the z 
coordinate as the direction of the waveguide axis and express the field 
vectors and the differential operator V as superpositions of transverse 
and longitudinal parts. The symbol t indicates the transverse direc- 
tions x and y. Thus, we have 

E = E, + B„ (7) 

H = H, + H £ , (8) 

and 

v = V, + e,~ (9) 

oz 

We use the notations e x . e„, and e* to indicate unit vectors in x, y, 
and 2 directions. 

The transverse part of the vector e-E is indicated by the notation 
e,E or, in component notation, 

e*E = e x (e xx E x + t xy E v + t xz E z ) (10) 

e„-E = e v (e vx E x + t vv E y + e vz E z ). (11) 

The longitudinal part is 

e z E = e z (e zx E x + UyE v + t zz E z ). (12) 

We may now separate Maxwell's equations into transverse and longi- 
tudinal parts. The transverse parts of (3) and (4) are 

V, XH 2 + e ; X^ = icoerE (13) 

oz 

V, X E, + e z X ^ = - 10,/xoH,. ( 14 ) 

oz 

Their longitudinal parts may be written as 

V t X H« - iw(«,-E« + 6„E 2 ) (15) 

and 

V t X Ei - - *o>moH z . (16) 
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The longitudinal parts of E and H follow immediately from (15) and 
(16), 

E 2 = J- V, X H, - — e z -E, (17) 

and 

H z = - J- V, X E,. (18) 

On the right-hand side of (17) and (18) appear only transverse com- 
ponents of E and H. It is important to distinguish between the single 
and double subscript notation of e. A double subscript, like e„, indi- 
cates a single tensor element of e, while a single subscript, like t z , is 
defined by (10) through (12). In particular, we have 

6,-E, = e z (e zx E x + e zv E y ). (19) 

We now use (17) and (18) to eliminate the z components of E and H 
from the transverse parts of Maxwell's equations (13) and (14), 

- J- V, X (V, X E,) + e 2 X ^ 
luino dZ 

= focrEi -— *.«,.£, + — tr(Vi X H,) (20) 

and 

V, X \^— V,XH,-- e,-E,l + e 2 X ^ = - iwmoH,. (21) 
\_lue zz e zz J dz 

These two vector equations represent four scalar equations. Once 
eqs. (20) and (21) are solved, the z components of E and H can be 
obtained by simple differentiation from (17) and (18). We have thus 
achieved a simplification of the original problem by reducing the num- 
ber of equations from six, in (3) and (4), to only four. 

The components of the e tensor are assumed to be functions of x, y, 
and z. The e tensor defines the wave-guiding structure. Because of the 
z dependence of e, eqs. (20) and (21) do not have mode solutions. A 
normal mode is defined as a solution of Maxwell's equations whose z 
dependence can be expressed by the simple function 

<r»". (22) 

Such solutions exist only if the dielectric tensor does not depend on 
the z coordinate. To construct solutions of the general eqs. (20) and 
(21), we consider solutions of simpler equations that are defined by 
a tensor e that is similar to e but is independent of z. The choice of e 
is obviously arbitrary and is determined by convenience. Using (22), 
we find from (20) and (21) the following equations for the normal 
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modes of the waveguide structure defined by 
- J— V, X (V, X £<{") - ijflj^e. X JC<?> 

IWflQ 



= foe, • £<?> - * - e, • 6, • £<f + J- Si • (V t X 3C<?>) (23) 



C;: e 2 



and 



V, X [ -A- V, X Xl? - -_- «,■£<? 

- #»e, X Sir = - iwMoWjp. (24) 

The subscript v indicates a mode label. Equations (23) and (24) admit 
an infinite number of solutions with different eigenvalues (propagation 
constants) 0i p) and different field vectors £< p) and 3Cj rt . Script letters 
indicate mode fields, while roman letters E and H are reserved for gen- 
eral field distributions. The modes are of two different types, guided 
modes whose fields are confined to the vicinity of the waveguide and 
radiation modes that extend to infinity in transverse direction to the 
guide. 23 Guided modes have discrete eigenvalues /Si"', while the eigen- 
values of radiation modes form a continuum. The superscript (p) 
stands for either ( + ) or ( — ), depending on the direction of wave 
propagation. A wave traveling in the positive z direction has positive 
(real) values /3£ +) , a wave traveling in the negative z direction has a 
negative (real) value /3£ -) . In isotropic media, we have the simple 
relations, 

#-> = - Pi + \ (25) 

£<,-> = 6<, +) , £<r' = - £L + ', (26) 

and 

XiT 1 = - XiP, 3Cir' = Xl+>. (27) 

General anisotropic media are more complicated, so that (26) to (27) 
do not apply. Modes traveling in one direction may be different from 
modes traveling in the opposite direction. 

III. ORTHOGONALITY RELATIONS 

The modes of anisotropic dielectric waveguides are mutually 
orthogonal. 2 For the purpose of deriving orthogonality relations, it is 
simpler to use Maxwell's equations in the form (3) and (4) instead of 
the form (23) and (24). Separating the z derivatives from the V 
operator, we write (3) for a mode labeled v and (4) for a mode labeled n, 

V, X Kl p) - #<»e. X 3C< P > = foi-8« (28) 

and 

V, X B« - #»e, X £<<» = - icouM™. (29) 
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Next, we take the complex conjugate of (28), multiply the resulting 
equation by 8™, multiply (29) by -K™*, add the two equations, and 
integrate over the infinite cross section (e is assumed real) : 

f [{&P'Vt X 3Qi p) * - 3C< P) *-V, X £<<" 

+ #«*ej«-e. X 3C< P) * + t'/S^aC^-e. X Efp}dxdy 

= - u* J jlSF-e-W - noKl^-X^dxdy. (30) 

The first two terms on the left-hand side of (30) can be expressed as 

- f f V, • (£<«> X Ki p) *)dxdy = - j (£<<" X 3C< P> *) nds. (31) 

The two-dimensional divergence theorem was used to convert the 
integral over the infinite cross section in the x-y plane to an integral 
over the infinite circle with outward normal direction n and line ele- 
ment ds. The integral on the right-hand side vanishes if at least one 
of the two modes is a guided mode. If both modes v and n are radia- 
tion modes, the integral vanishes in the sense of a delta function of 
nonzero argument. 2 Using this fact and a well-known vector identity, 
we can express (30) as 

W - 0£ p> *)//e 2 - (6j e) X Kf^dxdy 

= - co /" jlZ^-i-Z'f - noK^-Ki^dxdy. (32) 

Because of the symmetry of the e tensor, the following relation holds : 
£</»*.?.£<<» = g(«.e •£;<"*. (33) 

We take the complex conjugate of (32), interchange the superscripts p 
and q as well as the subscripts v and y. and, using (33), subtract the 
new expression from (32) with the result : 

W - 0?H / / e *-[ £ M ?) X 3C< P) * + £i'"* X K^dxdy = 0. (34) 

Equation (34) is the desired orthogonality relation. It is obvious that 
this expression holds also for isotropic media. However, in the iso- 
tropic case it is possible to use (25) through (27) to prove that each 
term in (34) must vanish separately. 2 For the general anisotropic case, 
(34) cannot be simplified further. We infer from (34) that the integral 
vanishes if 0® - /3£ p) * ^ 0. This means that the integral vanishes even 
in the case v = p if p and q indicate opposite signs, and a wave is 
orthogonal to its backward traveling counterpart (if /3i 8) is real) if 
orthogonality means vanishing of the integral in (34). 
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The integral in (34) expresses the total power flow if v = n and 
p = q. We may therefore use the orthonormality relation, 



//e*-[S^> X 3fC<?>* + £<?* X Knixdy 






=- 2s * a) -nw pa -« (35) 



to express mode orthogonality and normalization. The subscripts t 
indicating the transverse parts of the modes were added since the z 
components of the fields do not contribute to (35). P is a normalizing 
factor common to all modes that is used to adjust the arbitrary ampli- 
tudes of the normal modes. For real values of /3< 9) , we have 

s<<° = 1. (36) 

In this case, the sign of the integral is expressed correctly by the fact 
that jSj 8) reverses its sign if q goes from (+) to ( — ). For opposite signs 
of p and q, the right-hand side of (35) vanishes as required by (34) 
if @l 9) is real. For imaginary j3< 9> , (35) vanishes for q = p. The orthog- 
onality relation also holds for imaginary values of 0i"\ Imaginary 
values of the propagation constants occurs only for evanescent "radia- 
tion" modes. 2 - 3 In the case of imaginary /3£", the sign of the right-hand 
side of (35) is not certain. For this reason, we have introduced the factor 
s£ e) that must be adjusted so that P is a positive real quantity. This 
means that s^ 9) may have to be negative, s^ ff) = — 1. However, this 
case can arise only in connection with evanescent "radiation" modes. 
The 5„„ symbol in (35) indicates Kronecker's delta if both modes are 
guided. When one mode is guided while the other is a radiation mode., 
we have 6„ M = 0. If both modes are radiation modes, 5,„ must be 
interpreted as the Dirac delta function. 

IV. DERIVATION OF COUPLED-WAVE EQUATIONS 

Any arbitrary field distribution compatible with Maxwell's equations 
can be expressed as the superposition of all the modes of the idealized 
structure defined by the dielectric tensor e. Because the complete set 
of modes consists of a finite number of guided modes plus a continuum 
of radiation modes, we express the transverse parts of a general field 
by the expansion 



and 



E, - £ <*<">£<?> + £ /"V*>(p)Sp>(p)dp (37) 

v,p p JO 

H, = £ ai»Xl? + £ f X a^(p)3C^(p)dp. (38) 

V ,V p Jo 
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The longitudinal parts follow from (17) and (18). The superscripts 
assume the values ( + ) and ( — ) indicating waves traveling in positive 
and negative z direction. The first terms in (37) and (38) represent the 
contribution of the finite number of guided modes labeled v. The 
second terms indicated combinations of sums and integrals. The 
integration ranges over the entire region of continuous-mode labels p 
and includes radiation modes with real as well as imaginary values of 
/3 (p) (p). The summation symbol in front of the integral sign indicates 
that, in addition to modes traveling in positive and negative z direc- 
tion, various types of radiation modes exist and must be added to 
obtain the complete set of modes. For the purpose of deriving coupled- 
wave equations, the notation of (37) and (38) is too cumbersome. We 
use an abbreviated notation by omitting the integration sign, leaving 
it understood that the summation symbol includes summation over 
guided modes and summation as well as integration over radiation 
modes. We thus write 

E,= Za™S& (39) 

and 

H, = £ a«3C#>. (40) 

£$° and 3CS? are independent of z, but al p) is a function of z. Substitu- 
tion of (39) and (40) into (20) and (21) and use of the mode eqs. 
(23) and (24) leads to 

E(^ + ^a< p, )(e 2 X3C<?>) 

= E ai p) ■ iu(e t - h) ■ £<?' - i<» ( — W ** ~ ~ if u ) ' S#> 

v,p \ tzz *zz / 

+ (h et ~ hz 1 )' {Vt x ^l (41) 

and 

E(^ + WJ(e 2 X£!f) 

= - I^lv, X \±( — - ±-) V t X3C<?> 

_(.I e ,_J.A 8if) l 



(42) 



We take the scalar product of (41) with - £<?>* and of (42) with JC#>*, 
then we add the two equations, integrate over the infinite cross section, 
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and use the orthogonality relation (35). The result of this procedure is 
aid) _l_ flW* / ( i n w \ 

" i/sri U " " / 

= £ai P) f [ \-^^V*-(e t - e,)-£<f> 

v,p J J I 



-(=-£)■*]] 



dzdy. (43) 



On the left-hand side of (43), we have used a superscript (r). This 
notation is necessary to distinguish between the case of real and im- 
aginary propagation constants /3^\ If (3^ is real, we have r = q. If 
j8jS° is imaginary, as it is for evanescent radiation modes, we must choose 
for r the sign opposite to q. We now write (43) in the abbreviated form 

^- = - W + £ *# rt a< rt . (44) 

The coupling coefficient is defined by (43). We may eliminate the 
transverse magnetic-mode field vector from the coupling coefficient 
by using (17) (applied to the mode field) and the identity (which is 
obtained by partial integration), 

f f*y- ( v < X *)<*"*» = // (v < X ^P*)'Vdxdy. (45) 
The coupling coefficient can be expressed as 






K ,, )= i"\W\ f f \ S if.\(tiLl 

-^<£«- «.)•(£■«* + *) 

S'? 1 * + i„^D • [( l f - i)( ~ «* + sg } ) 

- (^-^\sl?~\\dxdy. (46) 



+ (?. 



V. IMPORTANT SPECIAL CASES 

In its complete form, (46), the coupling coefficient is very compli- 
cated. For isotropic media, where the e tensor degenerates to a multiple 

COUPLED-MODE THEORY 993 



of the unit tensor, (46) simplifies to the exact form, 

K ' ; " = 400^P i I (t ~ l) [ &r ■ £il " + t S "'' S "] dzdVl (47) 

in complete agreement with the well-known result. 12 

In many practical applications, the anisotropy of the dielectric 
medium is only slight, and the difference between the actual dielectric 
tensor e and the ideal tensor e is small. In that case, (46) can be sub- 
stantially simplified. A reasonable approximation of (46) for slight 
anisotropy and small values of e — I is 

K '-" = ii^P? / / £ "*' {t ~ ~ l) ■ tl " dxdy - (48) 

Note that the whole vectors of the electric mode fields enter (48) and 
not just the transverse or longitudinal parts. The approximation (48) 
is obtained by considering off-diagonal elements and differences be- 
tween diagonal elements of e and e as quantities that are small of first 
order. Products of two first-order quantities have been neglected. For 
readers who did not follow the detailed derivation, we repeat here 
briefly the definitions of symbols appearing in (46) through (48). The 
symbol w is the angular frequency of the electromagnetic field. The 
script symbols &i p) indicate the electric-field vectors of normal modes 
of an idealized waveguide that is defined by the dielectric tensor 
i = e(x, y), the subscript visa mode label, and the superscript (p) 
stands for either (+) or ( — ), indicating the direction of wave propaga- 
tion. The propagation constants p™ of the modes are labeled in the 
same way as the field vectors. The superscript r is usually identical 
with the superscript q. Only in the case of imaginary (3jp (this case 
happens for coupling to a nonpropagating radiation mode and is of 
little practical interest) does r indicate the sign opposite to q. Likewise, 
ajp = 1 for most cases of interest. Only for imaginary values of 0£ ] may 
it become necessary to choose s^ r) = — 1 to keep the power normaliza- 
tion coefficient P positive in (35). The asterisk indicates complex con- 
jugation. Subscripts t and z occurring in (46) and (47) refer to the 
transverse and longitudinal parts of the vectors to which they are 
attached. Similar subscripts attached to e are defined by (10) through 
(12) and (19). The dielectric tensor e = e(x, y, z) defines the actual 
waveguide (in contrast to the ideal guide that is only a mathematical 
fiction). The integrals are extended over the infinite transverse cross 
section of the guide. Equation (48) assumes the same limit as (47) if 
the dielectric tensor degenerates into a multiple of the unit tensor 
since, in the spirit of the approximation (48), we must use e zz /e tz = 1. 
Equation (48), even though it is only an approximation, is likely to be 
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of most importance in practical applications because of its simple form. 
For many practical problems, the approximation is justified and leads 
to sufficiently accurate results. 

VI. CONCLUSIONS 

We have derived coupled-wave equations representing exact solu- 
tions of the electromagnetic field problem for dielectric waveguides 
that consist of anisotropic materials whose dielectric tensor is a function 
of the 2 coordinate. The field of the general waveguide is expressed in 
terms of ideal modes of a hypothetical dielectric waveguide defined by 
a dielectric tensor whose elements are independent of the z coordinate. 
The main result of this paper is the expression (46) for the coupling 
coefficients. For many practical applications, the exact coupling coeffi- 
cient- can be approximated in the simple form (48) . 

The coupled-mode theory for anisotropic dielectric waveguides is 
essential for the solution of problems of mode propagation in integrated- 
optics guides with random or systematic irregularities. A particularly 
important area of applications are guides that are made anisotropic 
by an externally applied dc voltage or whose anisotropy is changed by 
such a voltage. Instead of an applied voltage, an acoustical wave may 
cause an anisotropic change of the refractive index of a dielectric wave- 
guide. These cases cannot be handled by the simpler isotropic coupled- 
mode theory, but require the extension to anisotropic media presented 
here. 
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